Detection of pairing correlation in the two-dimensional Hubbard model 
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Quantum Monte Carlo method is used to re-examine superconductivity in the single-band Hub- 
bard model in two dimensions. Instead of the conventional pairing, we consider a 'correlated pairing', 
(cjfCjM ) with Ci a = Cj CT (l — m-tr), which is inferred from the t-J model, the strong-coupling limit of 
the Hubbard model. The pairing in the d-wave channel is found to possess both a divergence like 
1/T in the pairing susceptibility and a growth of the ground-state pairing correlation with sample 
size, indicating an off-diagonal long-range order near (but not exactly at) half-filling. 

PACS numbers: 74.20.-z, 71.10.Fd 
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Ever since the discovery of high-T c cuprate supercon- 
ductors, there has been a controversy about its micro- 
scopic mechanism. Following a seminal proposal by An- 
derson JjJ ascribing it to a purely electronic origin, ex- 
tensive numerical and analytical studies have been done 
to explore a possibility of superconductivity in the Hub- 
bard model in two dimensions (2D). Apart from possible 
relevance to the high-T c mechanism, it has been a long- 
standing problem to clarify whether the Hubbard model, 
a simplest possible model for correlated electrons, can 
indeed superconduct [||. 

A number of studies suggest an effective attraction be- 
tween electrons in the model. For example, two holes are 
found to be bound in some range of the on-site repul- 
sion U near half-filling from numerical calculations for 
finite systems || Quantum Monte Carlo (QMC) re- 
sults show that the pairing susceptibility is greater than 
the susceptibility calculated without the interaction ver- 
tex between quasi- particles near half-filling |j| , indicating 
that the quasi-particles interact attractively. 

However, it is not known whether these attractive 
interactions do indeed imply a superconducting off- 
diagonal long-range order (ODLRO). We can argue 
against the ODLRO by referring to the fact that the 
susceptibility itself for finite U is suppressed from that 
for U = 0. In addition, the susceptibility diverges 
at most logarithmically with decreasing temperature T, 
which is exactly the T-dependence for the non-interacting 
case (inset of Fig. 0(a)). On top of these, the projector 
Monte Carlo (PMC) calculation for the ground state ||] 
does not show any sign of a growth with the system size 
in the equal-time superconducting correlation, indicating 
the absence of such an order at least for U < At, where t 
is the hopping integral between nearest-neighbor sites. 

Nevertheless, we can still raise two possibilities that 
the Hubbard model can indeed exhibit superconductiv- 
ity that has eluded detection so far. One possibility is 
related to the view that the pairing should be composed 
of 'quasi-particles' rather than bare particles. Although 



this is the point proposed by White et al. g , they have in 
fact calculated the conventional BCS pairing correlation 
(eifev j_ct, |C^) with c\ a the spin-a electron creation op- 
erator at site i by assuming implicitly that the bare elec- 
trons injected to the system will promptly turn into ap- 
propriately 'dressed' ones. However, as stressed recently 
by one of the present authors |J, the dressing may not 
take place so quickly and completely in a microscopically 
small system as usually encountered in numerical calcu- 
lations. In such a situation, we should dress the electrons 
'by hand' from the outset to detect superconductivity. 

The other, more exotic, possibility is that correlated 
systems such as the Hubbard model may possess, in the 
thermodynamic limit, a new type of ODLRO that can- 
not be detected by (cjf t^ct,, cL). In this case, we should 
characterize the condensate in terms of a certain opera- 
tor O distinct from CifCVj with an ODLRO in the sense 
that (O) ^ breaks the gauge symmetry. 

In either possibility, a search for the pertinent correla- 
tion function deserves investigation. We propose in the 
present Letter to look at (c^Cii jej, ^Cj^) as a candidate 

with c la =Cur(l — ni-a-) and n lrT = c\ a c la . In the strong- 
coupling limit, the Hubbard model can be reduced into 
the t-J model in which the electrons move with excluded 
double occupancies as represented explicitly by the oper- 
ator Ci„. Thus, in this limit, the conventional pairing cor- 
relation automatically turns precisely into {fy\fy> fy-^) 
Our proposal amounts that even for the Hubbard 
model with a moderate U, this can be a proper quantity, 
for which we have implemented QMC to obtain the cor- 
relation functions for the 'correlated pairing'. The results 
indeed turn out to exhibit an evidence for superconduc- 
tivity in the 2D single-band Hubbard model. 
The Hamiltonian TL for the model is written as 

n = -t 5^ (4 a a^ + h.c.) + uJ2 n iT n n- CO 

Through a large-scale QMC study, we have evaluated 
both the equal-time correlated pairing correlation in real 
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space for the ground state, 

p a (i,j) = (At(i)A a (i) + A a (i)At(j)}, (2) 
and the k — static correlated pairing susceptibility, 

^ = ^E / /S <A a (i ) 0)At(j,r))dr, (3) 
iV id ^ 

as a function of /3 = 1/T in an TV-site system, where 

= X! /c«( 5 )(ciT £ i+<U - c. a c 4+(5iT ). (4) 

As for the pairing symmetries we consider d x i_ y i (a = d) 
with f d (x, y) = Sy fi (6 x< i + 6 x ,-i) - 8 Xl0 (S y>1 + 8 y -x) and 
extended s (a — s) with f s (x,y) = 5 y $(8 x .i + 8 Xt -i) + 
8 x ,o(8 y ,i + Sy-l). 

The details of the QMC algorithm arc the follow- 
ing. We have used both the finite-temperature and the 
ground-state PMC formalisms. In both cases, we have 
employed the discrete Hubbard-Stratonovich transforma- 
tion introduced by Hirsch M . In the Trotter decomposi- 
tion, the imaginary time increment [At = r/ (number of 
Trotter slices)] is taken to be cither < 0.03 for finite- 
temperature susceptibilities or < 0.07 in the ground- 
state pairing correlation. We have adopted the stabi- 
lization algorithm used by several authors to investigate 
ground-state and low- temperature properties JTcj ] . Since 
the negative-sign problem prevents us from obtaining re- 
liable results for large interactions, we have increased U 
and/or decreased T up to the point where the ratio of 
the total sign to the total number of samples degrades to 
- 0.3. 

We first show the results for the finite-temperature d- 
wave pairing susceptibility, Sd, in Fig. |l|(a). Calculation 
is performed for an 8 x 8 lattice, where the band filling 
n is kept to 0.74 by adjusting the chemical potential for 
each temperature. For U — and U = 1 (in units where 
t = 1), Sd increases only logarithmically with (3(= 1/T). 
The logarithmic behavior is similar to that for the con- 
ventional susceptibility Sd shown in the inset. However, 
when we increase the repulsion to U = 3, Sd starts to 
increase more rapidly and can no longer be fitted to a 
logarithmic form within the error bars. Strikingly, the 
behavior fits rather well to a linear (oc 1/T) dependence 
for (3 > 5. The 1/T behavior is indicative of a long-range 
order, as exemplified by the 1/T divergence of the stag- 
gered magnetic susceptibility of the single-band Hubbard 
model at half-filling |Jll|. By contrast, S s in Fig. 0(b) 
is strongly suppressed when U is switched on, and barely 
increases with (3. 

When the band filling is decreased to n ~ 0.6, Sd 
shows much weaker dependence on either f3 or U (Fig. 
^|(a)). There the f3 dependence is so small that we can- 
not even decide whether it is linear or logarithmic within 
the error bars. If we increase the band filling to the 



half- filling (n = 1), on the other hand, S d is suppressed 
when U is switched on, where the susceptibility seems to 
assume a logarithmic /3-dependence as indicated in Fig. 
||(b). These results suggest that the strong divergence 
of Sd with [3 occurs only in a limited region of n near 
half-filling. 

Now we turn to the equal-time pairing correlation, 
p a (i,j) (Eq.(||)), in the ground state. We represent its 
real-space behavior by the quantity P a (R) , defined by the 
sum of p a (r, r + Ar) on a square at a distance \Ar x \ — R 
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R. This sum is intended to reduce statis- 



tical errors. In Fig. ||, we show the ratio of P a (R) for 
U = 3 to that for U = for both d x 2_ y 2 and extended s 
channels for a 10 x 10 system with n = 0.82. For the d- 
wave (Pd(R)), the pairing correlation at large distances is 
indeed enhanced by U, while the correlation is strongly 
suppressed for the extended s. There we should note 
that the results at R — half the sample size suffer from 
an effect of periodic boundary condition. 

An increase of the correlation function with U does 
not necessarily signify the occurrence of ODLRO as 
stressed by various authors. An appropriate quantity 
for this purpose is the k = Fourier component, i.e., 
J2jPa(0,j) — Y1rPc>(R) of the real-space correlation, 
and its system-size dependence in particular. This quan- 
tity should grow with the system size when a long-range 
order sets in. 

The QMC calculation is done for this quantity in 
two different sizes: a 10 x 10 lattice with 74 electrons 
(n = 0.740), and a 14 x 14 lattice with 146 electrons 
(n = 0.745). The choice of these sizes is made due 
partly to the fact that the band fillings are almost the 
same, but we also have paid special attention to the 
fillings satisfying the closed-shell condition, i.e., band 
fillings for which the non-interacting Fermi sea is non- 
degenerate. (The pairing correlation at open-shell fill- 
ings is inherently smaller, so that a mixture of results for 
close- and open-shell fillings would obscure the conclusion 
|]l2|.) The closed-shell condition also reduces greatly the 
difficulty arising from the negative-sign problem in the 
PMC calculation. The above two fillings are the closest 
closed-shell ones within the system sizes tractable in the 
QMC. Even for such close fillings, Pd(R) at very short 
distances (R = or 1) has a considerable filling depen- 
dence due primarily to the filling dependence of density 
(diagonal) correlations contained inherently in the pair- 
ing correlation. In order to focus on the longer-range 
(off-diagonal) part, we have evaluated P d by the sum of 
Pd(R) for R > 2. The values for P^'s thus obtained for 
the two system sizes almost coincide (with the difference 
being less than 0.0008) for U = 13 . The same is true 
as long as U is smaller than 1.5. 

When U is further increased, the long-range part of the 
correlation grows with U, where P d does concomitantly 
start to increase with the sample size as seen in Fig. H. 
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To be more precise, the amount of this increase exceeds 
the error bar for U larger than about 2. Although there 
still remains, to be rigorous, a possibility that this size- 
dependence contains an effect of the slight discrepancy in 
the filling, the result for the correlation may be taken as 
an indication toward the ODLRO, which reinforces the 
conclusion drawn from the T- dependence of the suscep- 
tibility. 

We expect that the long-range behavior of Pd and Sd 
continues to grow for larger U , in which the present corre- 
lated pairing correlation will eventually tend to the con- 
ventional one. We can then envisage a consistent pic- 
ture that encompasses a recent result [jL4j for the Mott- 
Hubbard regime of Emery's three-band Hubbard model 
p5[ , where an indication of a long-range order in the con- 
ventional d-wave pairing is detected for Ud comparable 
to the one-electron (anti-bonding d-p) bandwidth. 

As for the nature of superconductivity for moderate 
U, further investigation is necessary in the following two 
respects. First, we cannot at the present stage decide 
whether the correlated pairing either represents just a de- 
vice for detection of superconductivity in small systems 
or indicates a new ODLRO. Second, extensive studies for 
the t-J model in ID jl6| and 2D JllJ suggest that super- 
conductivity appears for a large enough superexchangc 
interaction J/t with apparently the same form of pair- 
ing as ours. However, these studies focus on a parameter 
region distinct from the present study. Remember that 
if the t-J model is regarded as an effective Hamiltonian 
for the Hubbard model for U — > oo, we end up with a 
perturbatively small J/t. Thus whether the mechanism 
for superconductivity in the t-J model with a large J/t is 
similar to the present case remains to be another future 
problem. 

In summary, we have shown from the results of the 
pairing susceptibility and the correlation function ob- 
tained in the QMC method that the 'correlated pairing' 
with a d x 2_ y 2 symmetry can become long-ranged in the 
2D Hubbard model with moderate U near half-filling. 

Numerical calculations were done on HITAC 
S3800/280 at the Computer Center of the University 
of Tokyo, and FACOM VPP 500/40 at the Supercom- 
puter Center, Institute for Solid State Physics, Univer- 
sity of Tokyo. For the former facility we thank Prof. Y. 
Kanada for a support in 'Project for Vectoralized Super 
Computing'. This work was also supported in part by 
Grant-in-Aids for Scientific Research from the Ministry 
of Education, Science, Sports and Culture of Japan. 



[1] P.W. Anderson, Science 235, 1196 (1987). 
[2] J.E. Hirsch, Phys. Rev. B 31, 4403 (1985). 
[3] J.R. Riera and A.P. Young, Phys. Rev. B 39, 9697 (1989). 
[4] A. Parola et ai, Physica C 162-164, 771 (1989). 
[5] E.Dagotto et at, Phys. Rev. B 41, 811, (1990). 
[6] SR. White et al, Phys. Rev. B 39, 839 (1989). 
[7] J.E.Hirsch and H.Q.Lin, Phys. Rev. B 37, 5070 (1988). 
[8] N. Furukawa and M. Imada, J. Phys. Soc. Jpn. 61, 3331 
(1992). 

[9] Y. Takada, J. Phys. Soc. Jpn. 65, 1544 (1996). 

[10] G. Sugiyamaand S.E. Koonin, Ann. Phys. 168, 1 (1986); 
S. Sorella et al, Int. J. Mod. Phys. B 1, 993 (1988); SR. 
White et al, Phys. Rev. B 40, 506 (1989); M. Imada and 
Y. Hatsugai, J. Phys. Soc. Jpn., 58, 3752 (1989). 

[11] Away from the half filling, a QMC result M shows that 
the SDW correlation is short-ranged. 

[12] This caution applies to previous QMC results as well. 

[13] If we include the R — 0, 1 components, Pd for two sizes 
would differ by 0.013 for U — 0, which would obscure the 
size dependence of the long-range part of the correlation. 

[14] K. Kuroki and H. Aoki, Phys. Rev. Lett. 76, 4400 (1996). 

[15] V.J. Emery, Phys. Rev. Lett. 58, 2794 (1987). 

[16] M. Ogata et at, Phys. Rev. Lett. 66, 2388 (1991). 

[17] E.Dagotto and J. Riera, Phys. Rev. Lett. 70, 682 (1993). 

FIG. 1. (a) The correlated d-wave pairing susceptibil- 
ity, Sd, for a fixed band filling n ~ 0.74 is plotted against 
P(= 1/T) for U = 3 (O) or U = 1 (A), while the solid curve 
represents the non-interacting case. The dashed curve is a 
logarithmic least-squares fit for U — 3 (for 4 < (3 < 6), the 
dash-dotted line a linear fit for (7 = 3 (for 5 < f3 < 10), and 
the dotted curve a logarithmic fit for U — 1 (4 < f3 < 10). 
The inset depicts the conventional d-wave pairing suscepti- 
bility, Sd, for comparison, where O represents the result for 
U — 3 with a logarithmic fit for 4 < j3 < 10, while the solid 
curve indicates the non-interacting case, (b) The correlated 
extended s-wave pairing susceptibility, S s , for U = 3 (O) 
with a logarithmic fit for 4 < j3 < 10 (dashed line), or for 
U = (solid line). 



FIG. 2. (a) Sd as a function of /3 for a fixed n ~ 0.6 with 
U = 3 (O) an< l the non- interacting case (solid curve), (b) A 
similar plot for n — 1 , where the dashed curve is a logarithmic 
fit for (7 = 3. 



FIG. 3. The ratio of the correlation function P a (R) for 
U = 3 to that for U — is plotted against R for the d (O) ° r 
extended-s (A) pairing. The system size is 10 x 10 with 82 
electrons (n = 0.82). 



FIG. 4. Pd is plotted against U. Number of electrons and 
system size are 146 electrons/14 x 14 (Oi n — 0.745) and 74 
electrons/10 x 10 (n,n= 0.740). 
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